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ABSTRACT 

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays 
(0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra 
from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. 
For the 1 1 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the 
luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation 
contributions to the soft X-ray luminosity (i.e., L aA gn and L x ,sv) for the remaining 24 Sy2s were estimated from 
the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a 
template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities 
derived from spectral fitting and the spread in the scaling factor, we estimated L Ai agn and L v ,sf from Monte Carlo 
simulations. These simulated luminosities agree with L v ,agn and T a ,sf derived from Chandra imaging analysis 
within a 3er confidence level. Using the infrared [Ne n] 12.8 // m and [O iv]26 ji m lines as a proxy of star formation 
and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total 
soft X-ray emission. This decomposition generally agrees with L x $f and L. v ,agn at the 3er level. In the absence 
of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star 
formation in galaxies hosting type 2 AGNs. 
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1. INTRODUCTION 

Supermassive black holes (SMBHs) and their parent galaxies 
co-evolve (Kormendy & Richstone 1995; Magorrian et al. 1998; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; Tremaine et al. 
2002; Haring & Rix 2004). In particular, observational and the- 
oretical work has established a link between accreting SMBHs 
(active galactic nuclei or AGNs) and host galaxy star forma- 
tion. Common mechanisms have been proposed for triggering 
star formation while fueling SMBH accretion, including galaxy 
mergers (e.g., Sanders et al. 1988; Hopkins et al. 2008) and sec- 
ular processes involving gravitational instabilities induced by 
spiral arms or galactic-sized bars (e.g., Kormendy & Kennicutt 
2004; Cisternas et al. 2011; Schawinski et al. 201 1). 

Disentangling star formation from AGN activity becomes a 
necessary endeavor to investigate SMBH and host galaxy co- 
evolution. This separation has been explored extensively and 
has included identifying optical and ultraviolet (UV) spectral 
signatures of starburst activity in AGNs (Cid Fernandes et al. 
2001), performing principal component analysis on infrared 
spectra of AGNs (e.g., Buchanan et al. 2006), and analyzing 
optical and infrared (IR) diagnostics that parameterize the 
relative contribution of AGNs to star formation, such as ratios of 
infrared and optical emission lines that indicate ionization field 
hardness (Genzel et al. 1998; Kewley et al. 2001; Kauffmann 
et al. 2003; Treyer et al. 2010; LaMassa et al. 2010, 2012) 
and strength of the polycyclic aromatic hydrocarbon features 
(Genzel et al. 1998; O’Dowd et al. 2009; LaMassa et al. 2010, 
2012). Here we extend the study of the interplay between AGN 
activity and star formation to the X-ray regime. 

Quiescent galaxies emit thermal and non-thermal X-ray 
emission. Hot gas energized by stellar winds and supernova 
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explosions radiates in soft X-rays (0.5-2 keV) whereas X-ray 
binaries (XRBs) and supernova remnants dominate the non- 
thermal emission above 2 keV (see Fabbiano 2006 for a review). 
Due to the relatively short life time of high mass X-ray binaries 
(HMXB s, r < 10 7 yr) and the short delay between a starburst 
and the generation of soft X-rays from hot gas, hard and soft 
X-ray emission can trace the instantaneous star formation rate 
(SFR). Indeed, X-ray emission from quiescent galaxies is well 
correlated with radio, infrared, optical, and ultraviolet SFR 
indicators (e.g., Ranalli et al. 2003; Rosa Gonzalez et al. 2009; 
Pereira-Santaella et al. 2011). Various X-ray SFRs have thus 
been calibrated in the literature based on the correlation of X-ray 
emission with far-infrared and radio luminosities (Ranalli et al. 
2003), X-ray emission from HMXBs and cumulative galactic 
X-ray point sources with far-infrared luminosities (Persic et al. 
2004; Persic & Rephaeli 2007), and X-ray luminosity with 
infrared plus ultraviolet emission (Pereira-Santaella et al. 201 1). 
However, low mass X-ray binaries, which trace stellar mass (M*) 
rather than star formation, also contribute to the X-ray emission 
and become the primary source of X-ray radiation in galaxies 
with low SFR/M* (Colbert et al. 2004; Lehmer et al. 2010). 

An active nucleus complicates using X-ray emission as an 
SFR indicator: accretion onto the SMBH dominates emission 
above 2 keV. Isolating the hard X-ray emission due to star for- 
mation then becomes prohibitively difficult when the nucleus 
is not resolved. However, soft X-ray emission can have com- 
parable contributions from scattered/reflected AGN emission 
and thermal emission from gas associated with starburst ac- 
tivity when the AGN is obscured or has very low luminosity. 
Determining the relative contribution of AGNs and star forma- 
tion to the soft X-ray emission can be achieved through high- 
resolution spectroscopy which resolves narrow emission lines. 
Various diagnostics can be used to distinguish between collision- 
ally and photoionized plasma, such as the width of the radiative 
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recombination continua (Liedahl & Paerels 1996) and ratios 
of forbidden, intercombination, and resonance lines in the 
Ovii triplet (Porquet & Dubau 2000). However, with current 
X-ray missions (i.e., Chandra and XMM-Newton) high- 
resolution spectra are only obtainable through grating observa- 
tions, necessitating long exposure times on X-ray-bright sources 
to achieve adequate signal-to-noise. How well can less time- 
intensive low-resolution CCD spectroscopy achieve the goal of 
disentangling starbursts from AGN activity in soft X-rays? 

Local obscured AGNs, Seyfert 2 galaxies (Sy2s), provide an 
ideal laboratory for answering this question. As these sources 
are relatively nearby and the accretion disk is hidden behind an 
obscuring “torus,” circumnuclear starbursts are visible. Hence, 
both AGN and host galaxy star formation can be modeled. 
Levenson et al. (2004, 2005) analyzed Chandra observations of 
two Sy2/starburst composites (NGC 5135 and NGC 7130) and 
demonstrated that a thermal component due to stellar processes 
was necessary to adequately model the soft X-ray emission. 
We expand this methodology to Sy2s in general to investigate 
the efficacy of modeling soft X-ray emission with a thermal 
component to describe host galaxy star formation. 

Our analysis focuses on two samples of Sy2s: an optically 
selected [O iii]5007 A sample (LaMassa et al. 2009) and a mid- 
infrared-selected 12 /xm sample (Spinoglio & Malkan 1989). 
Using Chandra and XMM-Newton observations, we decom- 
posed the soft X-ray emission into a star-forming and an AGN 
component. For Sy2s observed with Chandra that have signifi- 
cant host galaxy emission, we are able to resolve the point source 
associated with the AGN and can therefore remove it and analyze 
emission solely from the host galaxy, which we ascribe to star 
formation (L v ,sf)- For the Sy2s observed with XMM-Newton 
and the remaining Chandra sources, we fit the X-ray spectra 
with a combination of a thermal and power-law model. To ac- 
count for the presence of XRBs, we scale the soft thermal and 
power-law fluxes using star-forming galaxies from the XSINGS 
sample as a template (A. Ptak et al. 2012, in preparation), obtain- 
ing soft X-ray luminosity values associated with star formation 
(L v ,sf) and AGN activity ( L x A c,n )- We compare these estimates 
of L x sf and L x AC ,n with infrared (1R) spectral signatures that 
accurately describe star formation ([Nen]12.8 /xm luminosity, 
L [Ne ii] 1 Ho & Keto 2007 ; LaMassa et al. 20 1 2) and AGN activity 
( [ O iv]26 /im luminosity, L[o lv| ; Rigby et al. 2009; Diamond- 
Stanic et al. 2009; Melendez et al. 2008; LaMassa et al. 2010) 
to test the accuracy of this decomposition. 

2. SAMPLES AND DATA ANALYSIS 

The selection details for the [Om] and 1 2 /i m samples are 
presented in LaMassa et al. (2009) and Spinoglio & Malkan 
(1989) and are therefore only briefly mentioned here. The [O in] 
sample was selected from Sloan Digital Sky Survey (SDSS) 
Data Release 4 (Adelman-McCarthy et al. 2006). These sources 
live in the Seyfert 2 locus of the BPT diagram (Baldwin et al. 
1981) according to the Kewley et al. (2006) demarcation, and 
the sample is complete to a flux limit of 4 x 10 14 erg cm~ 2 s -1 , 
totaling 20 sources. Seventeen of these Sy2s were observed with 
XMM-Newton. The 12 fim sample is composed of Sy2s from the 
IRAS point-source catalog, is complete to a flux density limit of 
0.3 Jy at 12 /xm, and totals 31 galaxies. Twenty-eight of these 
sources have archival Chandra and/ or XMM-Newton data. 

In LaMassa et al. (2009, 2011) we discussed in detail the 
X-ray data reduction, using XAssist (Ptak & Griffiths 2003), 
and broadband X-ray (0.5-10 keV) spectral fitting of these 


sources with XSpec. Each galaxy was fit with an absorbed 
power law or when necessary, double absorbed power law, with 
the spectral indices initially tied together, simulating a partial 
covering geometry where a certain percentage of the transmitted 
X-ray emission is absorbed and the remainder scattered into the 
line of sight; Fe Ka emission, when present, was modeled with 
a Gaussian component. To accommodate the possible presence 
of starburst emission in soft X-rays (0.5-2 keV), we included 
a thermal model (APEC) in these spectral fits which describes 
emission from a collisionally ionized plasma. Here we revisit 
more detailed modeling of the soft emission. 

Since our aim is to study host galaxy star formation in 
tandem with AGN activity, we limit our analysis to the Sy2s 
where the addition of a thermal model component results in 
a detection of thermal emission, rather than an upper limit. 
As a result, NGC 291, CGCG 064-017, CGCG 218-007, and 
2MASX J 1 3463217+6423247 were dropped from the [Om] 
sample. We also exclude the three 12 /xm sources that were not 
detected in 2-10 keV X-rays: F08572+3915, NGC 5953, and 
NGC 7590. NGC 1068, part of the original 12 /xm sample, has 
soft X-ray emission dominated by AGN photoionization (e.g., 
Kinkhabwala et al. 2002) and has therefore been omitted from 
this analysis. In total, 38 Seyfert 2 galaxies remain: 14 from the 
[O m] sample and 24 from the 12 /xm sample. 

As noted in LaMassa et al. (2011), several 12 /xm Sy2s had 
both Chandra and XMM-Newton archival observations. We 
tested if differences in the spectra for an individual source 
between the two observatories were due to aperture effects, since 
Chandra has a higher spatial resolution than XMM-Newton, 
or variability. Our aim in that analysis was to constrain the 
properties of the AGN, so we focused on the Chandra data 
to isolate the point source, using XMM-Newton information 
to help constrain the spectral fit for non-variable sources. In 
this analysis, we favor larger aperture extraction areas as they 
encompass extended X-ray emission due to host galaxy star 
formation. For the Sy2s in the 12 /xm sample with both XMM- 
Newton and Chandra observations, the Chandra extraction 
area was adjusted to match the area from XMM-Newton or to 
encompass all visible extended X-ray emission. For the sources 
with just Chandra data, we increased the extraction area from 
the analysis presented in LaMassa et al. (2011) to ensure all 
extended X-ray emission is covered. In Tables 1 and 2, we list 
the aperture sizes used to extract spectra for the [O m] and 1 2 /xm 
Sy2s, respectively. 

The spectral fitting results for the [Om] and 12 /xm samples 
are summarized in Tables 3 and 4. The broadband X-ray spectra 
(0.5-10 keV) were fitted to most fully constrain the AGN 
contribution to the soft X-ray emission. Multiple spectra of 
the same source (from the three XMM-Newton detectors and/or 
multiple observations) were initially fit simultaneously with a 
multiplicative constant allowed to vary to account for differences 
in detector sensitivity. For variable sources, individual model 
parameters were fit independently. We imposed a lower column 
density (Ah ) limit of that due to our Galaxy along the line of sight 
to the source. In some cases, the best-fit absorption pegged at this 
lower limit and Ah was subsequently frozen to this value. For the 
sources better described by a double power-law model, which as 
noted above represents a partial covering geometry, the photon 
indices were tied together, with the normalizations left free 
to vary and an independent absorption attenuating the second 
power-law component (Ah, 2 ); the Sy2s best accommodated 
with a single power law can be identified by the blank entry 
in the A H ,2 column of Tables 3 and 4. The lower temperature 
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Table 1 

[O hi] Sample: Aperture Extraction Areas 3 


Galaxy 

z 

Aperture Radius (") 
PN/MOS1/MOS2 

Aperture Radius (kpc) b 
PN/MOS1/MOS2 

Mrk 0609 

0.034 

65/34/35 

44/23/24 

IC 0486 

0.027 

34/33/35 

19/18/19 

2MASX J08035923+2345201 

0.029 

17/20/18 

10/12/11 

2MASX J08244333+2959238 

0.025 

40/40/40 

21/21/21 

2MASX J101 8 1928+3722419 

0.049 

20/22/22 

20/22/22 

2MASX Jill 10693+0228477 c 

0.035 

18/-/- 

13 

CGCG 242-028 

0.026 

19/22/13 

10/12/7 

SBS 1133+572 

0.050 

20/16/17 

20/16/17 

Mrk 1457 

0.049 

19/13/13 

18/13/13 

2MASX J1 1570483+5249036 

0.036 

35/35/35 

26/26/26 

2MASX J12183945+4706275 d 

0.094 

-/30/31 

53/55 

2MASX J12384342+0927362 

0.083 

33/30/30 

52/47/47 

NGC 5695 

0.014 

31/22/17 

9/7/5 


Notes. 

a [O in] sources were observed with XMM-Newton only. PN, MOS1, and MOS2 refer to the three detectors on board XMM-Newton. 
b Conversion from arcseconds to kiloparsecs is based on the cosmology of Hq = 70 km s -1 M pc 1 , O vi = 0.27, and = 0.73. 
c Only PN spectrum fit since MOS1 and MOS2 spectra suffered from low signal-to-noise. 
d Source fell on chip gap in PN detector so only spectra from MOS1 and MOS2 were fit. 


limit ( kT ) was set at 0.1 keV. Ten Sy2 galaxies were well fit 
by a second APEC component, significant at greater than the 
3a level according to the /-test: 2MASX J08244333+2959238, 
NGC 424, NGC 1320, NGC 1386, NGC 4388, M-3-34-64, 
NGC 5194, Mrk 463, NGC 7582, and NGC 7674. This result 
indicates the presence of a two-temperature gas, similar to the 
results of X-ray observations of starburst galaxies (Dahlem 
et al. 1998; Strickland et al. 2004): one phase at kT between 
0.6-0.75 keV and a lower temperature gas at ~0.15 keV, with 
the exception of M-3-34-64 which has two hot gas components 
at 0.77+° 0 ° 0 2 2 keV and 2.85+° 0 3 3 4 4 . 

We list the observed soft fluxes in Tables 5 and 6: total 
0.5-2 keV flux (Fo. 5-2 kcv), power-law component of soft 
X-ray flux (F pow ), and APEC component of soft X-ray flux 
(Fa pec)- In the subsequent analysis, we omit the Sy2s with 
a poorly constrained APEC flux, i.e., the sources where the 
error is an order of magnitude or more higher or lower than 
the measured value (2MASX Jill 10693+0228477 and CGCG 
242-028), leaving us with 35 Sy2s. 

3. EFFECTS OF X-RAY BINARIES 

Before assigning “thermal” (Fapec) flux to starburst activity 
and power-law emission (F pow ) to AGNs, we need to account 
for the effects of XRBs which contribute to the non-thermal 
emission. As noted above, HMXBs indicate instantaneous star 
formation in the host galaxy and this effect is incorporated 
into X-ray SFR calibrations. Since XRBs boost the non- 
thermal flux, F pow will overestimate the AGN emission while 
Fapec underrepresents host galaxy star formation. To assess the 
contribution of XRBs to the observable emission, we model 
just the extended emission when the AGN is resolvable. For 
the remaining sources, we apply a correction to Fapec and F pow 
using results of X-ray analysis of star-forming galaxies from the 
XSINGS sample (A. Ptak et al. 2012, in preparation). 

Spatially resolved emission from the narrow-line region 
(NLR) can also potentially affect the results, overpredicting 
the luminosity attributable to host galaxy star formation. For a 
handful of 12 /tm sources, the major radius (F mal ) of the NLR 
was measured from Schmitt et al. (2003a), while we used the 


relations from Schmitt et al. (2003b) to estimate NLR F ma j 
values for the remaining 12 /tm sources and the [Om] sample. 
For the [O ill] Sy2s, the NLR radius ranges from 280 to 900 pc, 
with a median size of ~0.5 kpc. Similarly, the 1 2 /mi selected 
Sy2s have NLR sizes ranging from 80 to 1320 pc, with a median 
F ma j value of ~0.39 kpc. For all [O in] and a majority of 12 /im 
sources, the NLR is an order of magnitude or more smaller 
than the X-ray extraction radius. The NLR is thus likely not 
resolvable at these scales, so we only focus on quantifying 
the effects of XRBs on the observed emission. One notable 
exception is IC 5063, where the radius of the NLR (~1.3 kpc) 
comprises a significant fraction of the X-ray extraction radius 
(3 kpc). We will return to this Sy2 in Section 5. Another potential 
bias could be due to photoionized gas from the NLR being 
indistinguishable from collisionally ionized gas at low CCD 
resolution. However, as we show in Section 5, such an effect 
does not systematically affect our results. 

3.1. Chandra Imaging 

Due to Chandra's arcsecond resolution and the close prox- 
imity of the 12 /x m Sy2s, the central point source (i.e., AGN) is 
resolvable for a subset of our sample. To accurately constrain 
the AGN position, we ran wavedetect on 2-10 keV Chandra im- 
ages of the 12 /xm sample. Spectra were then extracted from the 
remaining emission, omitting the AGN identified by wavede- 
tect. Similar to above, the broadband X-ray emission was fitted 
with (an) absorbed power law(s) and an APEC component. We 
only include those sources where the spectra from the residual 
emission have a high enough signal-to-noise to be fit with the 
X 2 statistic (minimum grouping of 10 counts per bin), therefore 
dropping the sources where the X-ray emission is concentrated 
in the nucleus: NGC 424, TOLOLO 1238-364, F05189-2524, 
and NGC 5347. The summary of the spectral fitting for these 1 1 
sources is listed in Table 7, and the associated total soft X-ray 
flux from the unresolved emission (i.e., host galaxy emission 
after AGN removal) is listed in the second column of Table 6. 
To first order, the total soft X-ray flux from the unresolved emis- 
sion can be considered an estimate of host galaxy star formation, 
L a ,sf- La. agn is then estimated by taking the difference between 
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Table 2 

12 /xm Sample: Aperture Extraction Areas 3 


Galaxy 

z 

Observatory 

ObsID 

PN/MOS1/MOS2 

Aperture Radius (") 
PN/MOS1/MOS2 

Aperture Radius (kpc) b 

NGC 0424 

0.012 

XMM-Newton 

00029242301 

34/33/33 

8/7/7 



Chandra 

03146 

15 

3 

NGC 1144 

0.028 

XMM-Newton 

0312190401 

36/33/33 

20/18/18 

NGC 1320 

0.009 

XMM-Newton 

0405240201 

37/35/35 

6/6/6 

NGC 1386 

0.003 

XMM-Newton 

0140950201 

45/40/40 

2/2/2 



Chandra 

04076 

19 

1 

NGC 1667 

0.015 

XMM-Newton 

0200660401 

32/22/34 

10/7/10 

F05 189-2524 

0.043 

XMM-Newton 

0085640101 

41/31/33 

35/26/28 



Chandra 

02034 

40 

34 



Chandra 

03432 

40 

34 

NGC 3982 

0.004 

XMM-Newton 

0204651201 

34/37/44 

3/3/4 



Chandra 

04845 

34 

3 

NGC 4388 

0.008 

XMM-Newton 

0110930701 

36/35/43 

7/7/8 



XMM-Newton 

0110930301 

37/32/33 

7/6/6 



Chandra 

01619 

36 

7 

NGC 4501 

0.008 

XMM-Newton 

0112550801 

18/22/33 

3/4/6 



Chandra 

02922 

20 

4 

TOLOLO 1238-364 

0.011 

Chandra 

04844 

20 

5 

NGC 4968 

0.010 

XMM-Newton 

0002940101 

17/18/14 

4/4/3 



XMM-Newton 

0200660201 

21/18/18 

5/4/4 

M-3-34-64 

0.017 

XMM-Newton 

0206580101 

59/34/38 

21/12/14 

NGC 5135 

0.014 

Chandra 

02187 

20 

6 

NGC 5194 c 

0.002 

XMM-Newton 

0112840201 

65/42/49 

3/2/2 



XMM-Newton 

0303420101 

-/56/49 

72/2 



XMM-Newton 

0303420201 

-/61/46 

-/3/2 



Chandra 

00354 

40 

2 



Chandra 

01622 

40 

2 



Chandra 

03932 

40 

2 

NGC 5347 

0.008 

Chandra 

04867 

6 

1 

Mrk 463 

0.050 

XMM-Newton 

0094401201 

49/45/45 

49/45/45 



Chandra 

04913 

40 

40 

NGC 5506 d 

0.006 

XMM-Newton 

0013140101 

60/-/- 

9/-/- 



XMM-Newton 

0013140201 

54/-/- 

8/-/- 



XMM-Newton 

0201830201 

56/-/53 

8/78 



XMM-Newton 

0201830301 

52/-/64 

8/-/9 



XMM-Newton 

0201830401 

55/-/55 

8/78 



XMM-Newton 

0201830501 

68/-/98 

10/-/ 14 



XMM-Newton 

0554170201 

64/-/- 

9/7- 



XMM-Newton 

0554170101 

66/-/- 

10/-/- 

Arp 220 

0.018 

XMM-Newton 

0101640801 

20/21/25 

8/8/9 



XMM-Newton 

0101640901 

28/25/23 

11/9/9 



XMM-Newton 

0205510201 

34/28/22 

13/11/8 



Chandra 

00869 

40 

15 

NGC 6890 

0.008 

XMM-Newton 

0301151001 

17/13/15 

3/2/2 

IC 5063 

0.011 

Chandra 

07878 

12 

3 

NGC 7130 

0.016 

Chandra 

02188 

17 

5 

NGC 7172 

0.009 

XMM-Newton 

0147920601 

52/41/34 

8/7/5 



XMM-Newton 

0202860101 

61/59/60 

10/9/10 



XMM-Newton 

0414580101 

53/59/55 

8/9/9 

NGC 7582 

0.005 

XMM-Newton 

0112310201 

44/41/41 

4/4/4 



XMM-Newton 

0204610101 

60/43/47 

6/4/4 



Chandra 

00436 

43 

4 



Chandra 

02319 

43 

4 

NGC 7674 

0.029 

XMM-Newton 

0200660101 

19/18/16 

11/10/9 


Notes. 

3 For XMM-Newton observations, the extraction areas for the PN, MOS 1 , and MOS2 detectors are listed separately. 
b Conversion from arcseconds to kiloparsecs is based on the cosmology of Ho = 70 km s -1 Mpc -1 , = 0.27, and Oa = 0.73. 

c Source fell on chip gap on PN detector for ObsIDs 0303420101 and 0303420201; these spectra were not fit. 

d Dashes indicate that the spectra suffered from pile-up in that particular detector and were therefore not fit. See LaMassa et al. (201 1) for details. 
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Table 3 

[O in] Sample: APEC Model Parameters (Solar Abundance) 


Galaxy 

ATh .1 

( 10 22 cm " 2 ) 

m 

(keV) 

kT 2 

(keV) 

r 

Nn.2 

( 10 22 cm- 2 ) 

X 2 2 APECs 
DOF 

X 2 1 APEC 
DOF 

Mrk 0609 a 

0.04 

0 27 +005 
u - z/ -0.04 


i 77+0.05 
1 ■ 7 7 —0.04 



160 (203) 

IC 0486 

<0.06 

<0.16 


1 70 + 0.08 
— 0.07 

1 oo+O - 10 

1>uu _ 0.09 


636 (629) 

2MASX J08035923+2345201 b 

n fii+o.io 
uol -0.14 

< 0.12 


7 04 +I.O 5 
— 1.18 

46 - 7 -24. 6 l 


120 (95) 

2MASX J08244333+2959238 a 

0.03 

0 - 68 t° 0 0 o 6 

< 0.12 

| 54 +O .36 
l-J^_o.34 

15.52*2 4 

229(186) 

263 (188) 

2MASX J10181928+3722419 a 

0.01 

018-o°05 


2 63+ 0 - 64 
-0.69 



52.1 (61) 

2MASX Jllll 0693+0228477 c 

<0.62 

n 7 -2+0.15 
u - z ' 3 — 0.09 


2 04 +1-88 
z - U 4 - 0.88 



37.5 (36) 

CGCG 242-028 b 

r\ z;q+0.17 
0-69- 023 

<0.15 


0-31-°o 4 49 



87.0 (90) 

SBS 1133+572 

< 0.10 

0.8242°o 2 2 3 1 


3-08^i 

57 6 +45 - 4 
J 7 -30.2 


38.1 (48) 

Mrk 1457 

n /;/;+ 0 .i 2 
0 - 66 — 0.12 

0 14+0-03 
u - a ^-0.03 


1 29+ 137 
a - z *-1.14 

27.579 5 3 8 


35.2 (35) 

2MASX J1 1570483+5249036 

< 0.1 

0 17 + 0.03 
u ' 1 7 —0.05 


7 oc+0.23 
-0.34 

83 - 9 -27 1 9 


123 (76) 

2MASX J 1 21 83945+4706275 a 

0.02 

<0.24 


| Q5+0.70 
0.86 

87 - 2 -34 8 l 


15(19) 

2MASX J12384342+0927362 

<0.07 

<0.23 


7 17 +O. 3 I 
7 - 0.22 

29-372.6 


195 (164) 

NGC 5695 

<0.59 

0 ?4 +0 - 73 
u - z ^- 0.12 


7 ^c+ 0.68 

Z - J - 0.52 



75.7 (62) 


Notes. 

a Best-fit absorption same as Galactic absorption. This parameter was then frozen to the Galactic value. 
b Used c-stat. 

c PN only, MOS1 and MOS2 had low signal-to-noise spectra. 


L x s f and Uo. 5 - 2 kev for these Sy2s, with the associated error de- 
rived by scaling the error of total soft X-ray luminosity by the 
ratio of the AGN luminosity to the total luminosity. 

3.2. Star-forming Galaxy Template 

We have used the results of XSINGS to estimate the contri- 
bution of XRBs to the soft X-ray emission in host galaxies. The 
XSINGS project entails the analysis of 96 Chandra observations 
of 56 SINGS galaxies (Kennicutt et al. 2003). SINGS is designed 
to cover a wide range of galaxy properties such as morphologi- 
cal type, SFR, stellar mass, and metallicity, and the main aim of 
XSINGS is to study the relationships between X-ray emission 
(XRBs, nuclear activity, and hot interstellar medium) and these 
properties. In A. Ptak et al. (2012, in preparation) the spatially 
averaged X-ray emission of XSINGS galaxies is discussed and 
compared to the summed XRB luminosities for galaxies with 
significant detections (XRB catalogs for all XSINGS galaxies 
will be presented in Jenkins et al. 2013, in preparation). XRBs 
emit as a power law, but provide a measure of the instantaneous 
SFR along with thermal emission in soft X-rays. We therefore 
have 

i.t.SF = Fapec + Fxrb (1) 


Fa, AGN — Lpow — Lx RB ■ (2) 

We further assume that a certain fraction of X-ray emission from 
XRBs contributes to host galaxy star formation: 


Lx RB — R X L v,SF- 


(3) 


Using the 22 pure star-forming galaxies in XSINGS (i.e., non 
Sy2s and non LINERs), we have calculated the mean ratio of 
X-ray emission from resolved non-nuclear point sources (i.e., 
XRBs) to total X-ray emission in the 0.5-2 keV band (R), 
deriving R = 0.51 ± 0.26. Solving the above equations, we 
obtain 


L x , sf = 


L APHC 

1 - R 


(4) 


■^jc.AGN 


= L 


pow 


R 


Eapec 
\ — R 


(5) 


To account for uncertainties in the luminosities derived from 
spectral fitting (i.e., L A pec and L pow ) and the wide spread in R 
values, we have performed Monte Carlo simulations to derive 
L x sf and L x agn- We have drawn 1000 random values for 
Lapec, /-pow, and R from a Gaussian distribution centered on 
the best-fit values of Lapec and L pow (and the mean value of R), 
with a corresponding to the errors derived from spectral fitting 
(and standard deviation on R). Using these simulated parameters 
and Equations (4) and (5), we calculated 1000 simulated values 
for L x sf and L Ai agn- The mode of these distributions is then 
taken as L x sf and L x agn, with the 68% confidence interval as 
the associated errors. 

For the subset of Sy2s with high signal-to-noise Chandra 
spectra of the soft unresolved emission, we compare the results 
of L Xj sf (E.i.agn) derived via Chandra imaging analysis and 
the simulation approach described above. Here, L x sf Chandra 
(L x , agn, C handra) represents the soft X-ray luminosity due to star 
formation (AGN) based on spectral fitting of just the unresolved 
emission, while L XjS f, xsings (L x , agn, xsings) is the estimate of 
the star formation (AGN) component of the soft X-ray flux 
calculated via simulations. 

To account for errors in both parameters, we use a Bayesian 
approach to linear regression (Kelly 2007). Both here and in 
Section 5, the high and low error bars on the X-ray flux were 
averaged to provide symmetric errors for the linear regression 
routine. We plot L x , S f,xsings as a function of L XiSF , chandra 
and L Xi agn, xsings against L XjAG n, C handra in Figure 1. The gray 
shaded regions enclosed by the dashed lines delineate the 3er 
confidence level on the regression line, based on the median 
of the posterior distribution of the slope, intercept, and mean 
distribution of the independent variable, with the dot-dashed 
line denoting the line of equality. Though the T x ,sf, xsings values 
are somewhat systematically higher than i x ,sF,chandra, they do 
agree within the 3er contours. The agreement among these 
independent methods suggests that using XSINGS galaxies 
as a template to correct the Lapec and L pow values for the 
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Table 4 

APEC Model Parameters (Solar Abundance 3 ) 


Galaxy 

Wh,i 

m 

kT 2 

r 

A/ [ 1.2 

/ * 2 2 APECs 

X 2 * 1 APEC 


( 10 22 cm' 2 ) 

(keV) 

(keV) 


( 10 22 cm' 2 ) 

DOF 

DOF 

NGC 0424 

0 o4 +0 03 
u,u — 0.03 

0 69 +008 
U -O^_0.05 

<0.13 

a ni + 0.21 
Z-ZA — 0.28 

32.6^ 6 5 o 

246 (178) 

267(180) 

NGC 1144“ 

0.06 

o ^7+0.29 

U - 3 7 -0.06 


1 q 1+0.37 
A - yA -0.24 

47.(t 3 3 5 8 2 


175 (149) 

NGC 1320 

<0.08 

0 1 5+0.03 
U-A -0.03 

n 73+0-09 
u - ' J - 0.02 

9 Q7+0.27 
z - y/ — 0.21 

48 - 0 -15.5 

233 (188) 

279(190) 

NGC 1386 

o Q 4 + 0.02 
uu ^- 0.01 

0 67 +0 03 
U- ° '— 0.03 

0 13+0.02 
l-'-iJ-O.Ol 

9 77 +O. 12 
Z ' /Z — 0.14 

31 4 +22 - 9 
1 ‘*—11.5 

398 (332) 

435 (334) 

NGC 1667 b 

0.05 

n tio+ 0.07 
u - JJ -0.04 


9 1 q+ 0.34 
+-io_0.37 



49.8 (38) 

F05189— 2524 b 

0.02 

0.1 


2 05 +014 
Z -^_0.14 

6 77 + 0.44 
' ' -0.42 


501 (374) 

NGC 3982 b 

0.01 

0 29 +0 03 
u - zy -0.03 


9 39 + 0.18 
Z -J^_0.15 

40 3 +25 - 5 


174(160) 

NGC 4388 (XMM) b 

0.03 

0 60 +0 04 
U.OU-003 

0 15 +002 
U,A -0.03 

1 25 +012 
A ’ ZJ — 0.12 

24 3 +L1 
z ^" 5 -l .0 

510(495) 

580 (498) 

NGC 4388 (Chandra) b 

0.03 

0 60 +0 04 
0.04 

0 15+0.04 
U,A - 0.02 

O 98+0.29 
^ - 0.30 

25.6% 

210(166) 

269(168) 

NGC 450 l b 

0.03 

0 36 +0 05 
u - 3O -0.03 


1 C 9 +O. 3 O 
1 • JZ — 0.29 



94.2 (102) 

TOLOLO 1238-364“ 

<0.19 

0.32!_° 0 ^ 


2 45 + 0.33 
z -* J -0.36 



91.6(105) 

NGC 4968 b,c 

0.08 

°- 68 -d.i 2 


1 79 +O.I 8 
1 ’ 7 z — 0.21 



325 (268) 

M-3-34-64 b 

0.05 

7 q^+0.34 
^•O->_0.38 

n 77 + 0.02 
u - ' — 0.02 

3 24 +019 
— 0.23 

53.6 + _™ 

780 (492) 

848 (493) 

NGC 5 1 35 b 

0.05 

0-73To3 


2 34 + 0.07 
Z - J -0.08 



166(114) 

NGC 5194 

0 Q 4 + 0.01 

u - u ^_o.oi 

O-ISToi 

0 62 +0 - 01 
u.oz_ 00 i 

7 1 8+0. 14 
J-iO— 0.14 

l°- 2 -0.75 

1560(1291) 

2021 (1293) 

NGC 5347 b,c 

0.02 

< 0.21 


1 C 9 +O. 3 O 

1 dd_o.29 

32 6 +24a 
■ 3Z ‘ — 19.6 


69.9 (82) 

Mrk 463 b 

0.02 

0 74 + 0.04 
u - /i+ _ 0.02 

n 20 +0 06 
U-ZAA — 0.04 

1 76 + 0.13 
1 • /D — 0.16 

28 -6-3 5 4 

365 (316) 

392 (318) 

NGC 5506 d 

0 1 1 +001 
u - 11 - 0.01 

0 74 + 0.05 
u - '*-0.05 


171 + 0.02 
1 /A - 0.01 

2 £ 9 + 0.02 
z -° -0.03 


2646 (2380) 

NGC 5506 e 
NGC 5506 f 

f) 17 + 0.01 
U - A 7 - 0.01 

0 94+0.39 
0.24 


1 q 9 + 0.02 

l-OJ- 0.02 

9 on+ 0.01 

Z.5U-0 02 


3982 (3137) 

NGC 5506® 
NGC 5506 h 

0 12 +001 
u - AZ -o.oi 

n c'T+o.os 
-0.03 


1 7 A+ 0.01 

1 -/O -0 00 


» 

- 

NGC 5506' 

” 

0 96 +005 
u -^°_0.05 

” 

” 

” 

” 

” 

Arp 220“’ b 

0.04 

0-791° 0 Z 


0 95+0.26 
u - y -0.24 



309(302) 

NGC 6890“ 

< 0.22 

0-78 + _° 0 f g 


9 90 +O .88 
d-^o_074 

27 - 4 -i 8 i 4 3 


164(148) 

IC 5063 b 'i 

0.06 

o-eo^o. 1 ! 0 ! 


1 o^+0.16 

i.» 2_0 2 i 

2 1 .oiV . 2 


198(141) 

NGC 7130“ 

0 O 6 +0 03 
uuo - 0.02 

0.632°o°o 3 3 


2 43 + 0.25 
z - qo — 0.18 

75 4 + 55.6 
38.1 


334 (240) 

NGC 7172 b 

0.02 

0-28To4 


1 c£+ 0.02 

1 -50-0 03 

7 9 £+ 0.10 
7.3o_o 10 


2074(1746) 

NGC 7172 b 

” 

0 26 +0 02 
U ' ZO - 0.02 


1 54 +O.O 2 
A - J *-0.03 

8.133Y1 

” 

” 

NGC 7582 ( Chandraf 

0.01 

0 72 +0 05 
u - /z -0.05 


1 94 +O.H 

A - y *-o.io 

24.674 


355 (301) 

NGC 7582 (XMM) b 

0.01 

/a i O+0.02 
U - 1S -0.04 

0 72+0- 01 
u - /z -o.oi 

1 7 ^+ 0 . 05 
1 • '^—0.03 

27.°74 

1438 (886) 

1586 ( 888 ) 

NGC 7674 b ’“ 

0.04 

< 0.11 

n £7+0.07 
0- 6 / _o.o5 

0 62 +0-47 
U-OZ — 0.43 

89.074 ° 2 

66.6 (70) 

113(72) 


Notes. 

a All abundances frozen to solar except for Arp 220 which has an abundance of 0.17 *qqj. 

b Best-fit absorption same as Galactic absorption. This parameter was then frozen to the Galactic value. 

c Used c-stat. 

d ObsIDs 0201830201, 0201830301, and 0201830401. APEC and first power-law normalizations fit independently. 

e ObsID 0013140101. 

f ObsID 0201830501. 

8 ObsID 0013140201. 

h ObsID 0554170201. 

‘ ObsID 0554170101. 

1 Used pile-up model. 


presence of XRBs provides a reasonable estimate to decompose 
the soft X-ray emission in the absence of resolvable nuclear 
emission. 

In the subsequent analysis (i.e., Section 5), we use 

L x, SF, Chandra and f'x.AGN, Chandra US x , S h and L x AON fof the 1 1 

Chandra sources with high signal-to-noise. L a ,sf and L a ,agn 
for the remaining 24 Sy2s are estimated via the Monte Carlo 
simulations outlined above. These values are listed in Tables 8 
and 9 for the [Om] and 1 2 /tin samples, respectively. 


4. CONSTRAINTS FROM INFRARED DATA 

We compare our X-ray decomposition with information 
gleaned from high-resolution infrared spectroscopy, using di- 
agnostics that cleanly trace the separate processes of AGN 
and star formation activity. The flux of the [O i v ] 2 6 /ini line 
has been shown to accurately describe intrinsic AGN flux 
(Rigby et al. 2009; Diamond-Stanic et al. 2009; Melendez et al. 
2008; LaMassa et al. 2010) while the [Nen]12.8 /xm emission 
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Table 5 

[O in] Sample Fluxes (10 -14 erg cm -2 s _1 ) 


Galaxy 

Fo. 5— 2keV 

Fa pec 

^*pow 

Mrk 0609 

W-l-wa 

4 82 +214 
^• 6z -2.15 

74.32 2 2 9 9 2 2 

IC 0486 

35-0-O4 

O+S^oil 

34-6 + 4 7 8 s 3 

2MASX J08035923+2345201 

1 ?4 +0-86 
1 -Z — 0.71 

0.83!°o 8 7 4 o 

0 41 +O .17 
u -^ 1 -0.14 

2MASX J08244333+2959238 

r 25 + 1 - 26 

J - ZJ -1.58 

3 04+1-04 
— 1.47 

771 +0.71 
Z-ZA — 0.59 

2MASX J101 8 1928+37224 19 

7 97+0.92 
z - z/ -0.62 

1 04+0-80 
1 -U — 0.43 

1 77 +O .46 
A,z — 0.44 

2MASX Jllll 0693+0228477 

l 44+419 
1 — 0.59 

0-422 4 0 1 ] 6 

1 02 +1 ‘ 31 
1 -UZ — 0.46 

CGCG 242-028 

| 1-7+2.52 
1,1 —0.75 

O.80 + _ 2 0 5 7 °3 

0 217 +O .33 
U - J/ _0.19 

SBS 1133+572 

3 46+ 1 - 27 
0.82 

0-522°rii 

2 04 +I .22 
z ' y -0.74 

Mrk 1457 

2 on+10.3 
z,yu -1.82 

2.362 1 1 °79 

0 54 +O .66 
^•^_0.32 

2MASX J1 1570483+5249036 

-3 0 0 + 1.31 

0 O1+1-09 
u - yi -0.41 

2 Q2+0-72 
z -^ z — 0.51 

2MASX J121 83945+4706275 

1 54 + 2 - 15 
1 ' J -0.54 

0.55^j_3® 

0 oq+ 0.47 
u>yy -0.44 

2MASX J12384342+0927362 

6 23 +2-81 
°- z -0.76 

0 . 86 + _ 2 0 f 9 

5 37 +I.I 2 
J * J, -0.65 

NGC 5695 

2 C 7 + 5.37 
z - J/ -0.90 

0-64^ 

1 Q 7 + 3.20 
1 ■■'■^—0.71 


reliably tracks star formation in both quiescent and active sys- 
tems (LaMassa et al. 2012). We therefore use these fluxes to 
decompose the soft X-ray emission into a star-forming and an 
AGN component by solving Lo.5-2keV = a x L[Nen] + P x Zqoiv] 
using ordinary least-squares multiple linear regression on both 


the 12 /im and [Om] samples, obtaining a = 1.07 and ft = 
0.16 (see LaMassa et al. 2010, 2012 for calculated values of 
L[Oiv] and i[Nen] for these samples). Here Lo.5-2keV represents 
the total soft X-ray emission (i.e.. Column 1 in Tables 5 and 6). 

Figure 2 illustrates the results of this decomposition. We 
consider the X-ray variable Sy2s, NGC 5506, and NGC 7172 
as independent data points, but we do average the soft 
X-ray flux for the latter two NGC 5506 observations, which 
are consistent. Though we have only one IR measurement for 
each of these sources, we do not expect the infrared emission to 
vary as much as the X-ray emission. We color-code the sources 
in Figure 2 according to sample: cyan triangles represent the 
12 /rm Sy2s (with the X-ray variable sources connected by a ver- 
tical line), blue diamonds denote the subset of 12 /im sources 
with high signal-to-noise Chandra imaging of the unresolved 
emission, and red squares mark the [O m]-selected sources. The 
overplotted line indicates where the two quantities are equal. 
With the exception of a handful of outliers, a relatively good 
agreement is present between the observed soft X-ray emission 
and the IR decomposition. We note that this consistency holds 
for both the 1 2 /i m and [Om] samples, with no apparent sys- 
tematic offset introduced from different sample selection tech- 
niques. In the following analysis, we therefore use a x LjNen] 
and ft x L[Oi V ] as an independent estimate of the starburst and 
AGN contribution, respectively, to the soft X-ray emission. 


Table 6 

12/r.m Sample Fluxes (10 ~ 14 erg cm -2 s _1 ) 


Galaxy 

fb.5-2keV 

F D . 5—2 ke V , extended 3 

Tapec 

^pow 

Comments 

NGC 0424 

31.5 + _ 12 f 3 



20 9 +3 - 26 
zu>y -3.43 


NGC 1144 

8.39 41 , 25 ! 


2-87 4 ° 1 ' 3 9 

c C7+0.90 
J - JZ — 0.81 


NGC 1320 

28 4 +15 - 8 
zo -^— 5.46 


10.5 4l 4 5 92 

| 7 9+2-38 
1 />y -2.73 


NGC 1386 

7c 7+6.35 
ZJ ' Z — 3.99 

8.48 4l i 41 6 

i i 7+6.19 
1 a>/ -3.83 

1 3 5+1-43 
A * 5,J — 1.12 


NGC 1667 

7 31 +l45 
/ ' Ji — 1.52 


7 77+1.03 
■ 3 ’ zz — 1.10 

4 -09t 1 iSs 


F05 189-2524 

12 1 +2 - 02 
AZ ' A - 1.00 


7 77+I.88 
Z-ZZ — 0.68 

9.88 40 0 74 


NGC 3982 

9 14+O.9O 
v - 14 -0.96 

7.47^.2.49 

3-39 4 q 73 

tr 7^+0. 64 
J - /J -0.62 


NGC 4388 

28.5 41 7 6 q° 

23.3 46 6 6 4 

17.8 4l 6 5 6 4 

10 7+ 3 -15 
AU - 7 — 2.18 

Chandra and XMM-Newton observations averaged 

NGC 4501 

7-24 + . 1 1 5 3 7 o 

5 47+1.42 
d -^ / -1.47 

7 07+ 1.06 

-0.95 

3-41-0.88 


TOLOLO 1238-364 

14.2+f 


5 40+5.OI 
J-HO -3.80 



NGC 4968 

5-85L 1 1 1 o 4 5 


1 77+O.67 
1 -Z — 0.60 

4 50+0.93 
4>:)5 -0.86 


M-3-34-64 

50.8L 6 7 3 8 7 i 


00 4+4.72 

DJ.H-_7 29 

1 7 4 + 4 - 28 
1 ' ^-2.80 


NGC 5135 

34.5 + _ 1 1 “ 

16.4 474 

17 0 +1 ' 23 

i /.u_j 27 

17 5 +112 
1 '-^-l.lO 


NGC 5194 

69 6+6.82 

oy -°-6.07 

56.1 48 7 ° 9 

57 7+6.71 
JZ ‘ 7 —5.97 

16 9 +119 
ADy -l .09 


NGC 5347 

2 94+0-87 
Z- ^ — 0.46 


0 50+0.80 
U.JU — 030 

2.44 - 34 


Mrk 463 

1 4 0+2.68 
ih.7-272 

604 + _ 9 1 4 ° 

7 50 +2-59 
/-JU -2.62 

6.80 4 ° 0 ™ 


NGC 5506 

308 45 5 6 3 


6 07+ 4 - 48 
°’ U/ — 4.55 

302 45 5 j 

ObsIDs 0013140101, 0201830201, 0201830301, 0201830401 

NGC 5506 

444 4 j 7 7 


0 

444 + 37 

ObsIDs 0013140201 and 0201830501 

NGC 5506 

457 43 4 3 


10.5 + _‘if 2 

446 437 

ObsIDs 0554170201 and 0554170101 

Arp 220 

4.90 4l 1 1 3 4 

4 3R +111 
4 -5°-0.77 

3.05 4t >- 9 2 5 4 

| or +0.71 
l-OJ-0.51 


NGC 6890 

1 1 5 +619 
iA - D -3.76 


7 69+ 1 - 64 
z - oz _ 1 .40 

0 0 0+5.97 

O.OO — 3 49 


IC 5063 

22 9+3.12 
ZZ- — 3.18 

4 79+1-15 
^■ /y -1.20 

4 57 +1 - 58 
^• J/ -1.58 

1 0 7+2.69 
18. -5-2 76 


NGC 7130 

23 6 +2 - 63 
z ’ 3 “°— 1.89 

o cc+3.65 

°.jj _1 20 

10 3 +1 - 24 
AU -0.96 

13.3 423 6 2 


NGC 7172 

20 3 +2 - 28 
zu -- 3 _ 2.45 


2 05 +0 - 50 
Z ' UJ -0.51 

1 0 7 + 2.22 

1 o-5_2 40 

ObsID 0414580101 

NGC 7172 

1 1 o+u.75 

1 l-O— 0.72 


2 06 +0-27 
z,u — 0.28 

9 74+O.7O 
y>/ ^-0.67 

ObsIDs 0147920601 and 0202860101 

NGC 7582 

40.3 4 3 2 

30.8 43 2 2 

1 5 2 +2-5 
A * A-Z — 3.1 

25 l +2- ° 
zj - a - 2.0 

Chandra and XMM-Newton observations averaged 

NGC 7674 

17 5 +2 - 91 
1 ' -4.95 


13 6 +217 

A -3. 0-4 74 

3 90+1.94 
Jmyu -1.42 



Note. a Sy2s observed with Chandra where the unresolved emission after removal of the AGN had high enough signal-to-noise for adequate spectral fitting. 
fj).5-2keV, extended flux corresponds to this extended emission from the host galaxy (i.e., with the AGN removed). 
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(a) (b) 


Figure 1. Comparison of using Chandra observations to decompose the soft X-ray emission into star formation and AGN components (L x sf , Chandra and L x agn, Chandra, 
respectively) with estimates of starburst and AGN activity calculated via simulations using XSINGS data of normal galaxies and luminosities derived from spectral 
fitting (i.e., Ac.sr. XSINGS, Ta,agn, XSINGS )- The gray shaded regions enclosed by dashed lines indicate the 3 a confidence interval from a Bayesian linear regression fit 
(Kelly 2007). The dot-dashed lines denote where the two quantities are equal. These two methods of estimating the star formation and AGN contributions to the soft 
X-ray emission are consistent. 

(A color version of this figure is available in the online journal.) 


Table 7 

Chandra Spectral Fits to Unresolved Emission Only 2 


Galaxy 

1V h ,i 

( 10 22 cm” 2 ) 

kTi 

(keV) 

(keV) 

r 

Mi , 2 

( 10 22 cm” 2 ) 

X 2 2 APECs 
DOF 

X 2 1 APEC 
DOF 

NGC 1386 

<0.05 

0-58!°dT 3 


t zio+0.38 

48 9 +32 - 0 
12.4 


79.0 (54) 

NGC 3982 b 

0.01 

0 - 20 + _ 0 o°o 6 6 


2-04 + _° 0 f 8 



46.5 (25) 

NGC 4388 b 

0.03 

0 60 +0 04 
uou -0.04 

0 ifi+0.04 
U-iO -0.02 

°-98!°o 5 5o 


85.4 (83) 

123 (84) 

NGC 450 l b 

0.03 

0-37 + _°o 2 o 2 7 


2 26 +0 - 44 
^• zo _0.50 



31.2(38) 

NGC 5135 

<0.13 

0-84^® 

0 3Q + 011 
u - z, ^-0.09 

2 94+0.33 
0.31 


62.8 (60) 

77.9 (62) 

NGC 5194 b 

0.02 

0 6 O +0 01 
u - ou - 0.01 

fi i O+0.03 

U.Io_oo 2 

T 1 o+0. 10 

J,lD —0.09 

7 92 +1 - 38 
' '^-1.08 

525 (352) 

616 (353) 

Mrk 463 

<0.05 

0 62 +0 05 
uoz - 0.07 


2 50+0-75 
Z.OU-04? 

34.6 * 1 7 ° 7 5 


86 (71) 

Arp 220 b 

<0.04 

0 73 +0 - 07 
U-/ - 0.12 


2 1 q+ 0.44 
z ’ lu -0.36 



69.3 (57) 

IC 5063 b 

0.06 

0 34 + 0 09 
u.^_ 0 .04 


1 51 +0 - 55 
1 • J — 0.59 

28.5!_ 6 5 2 o 


71.5 (57) 

NGC 7130 

0 05 +0 07 
u.ud_oo3 

0 47 + 0.07 

u **'-o.io 


2 24 +0,37 
z ’ z ^-0.29 



58.2 (50) 

NGC 7582 

0 03 +0 02 
U-U - 0.02 

0 73 + 0.04 
U-/ -0.05 


9 97+0-23 
' -0.20 

41.82 1 8 °2 


142(128) 


Notes. 

2 All abundances frozen to solar. 

b Best-fit absorption same as Galactic absorption. This parameter was then frozen to the Galactic value. 


Table 8 5 . RESULTS : COMPARISON OF SOFT X-RAY 

[O hi] Sample Star Formation and AGN Luminosities DECOMPOSITION WITH IR 


Galaxy 

log(L x ,sF) 

l°g(L x ,AGN) 

Mrk 0609 

41.3 8L'o39 

42.28*0Q 6 6 

IC 0486 

40.39^042 

41.83*0^4 

2MASX J08035923+2345201 

40.45+° 0 8 3 7 3 

39. 64*0 3 3g 

2MASX J08244333+2959238 

40.8U° 0 S 3 0 3 

40. 1 8* 0 3 3 q 

2MASX J 101 8 1928+3722419 

40.99!°o 7 3 2 

40.54* 0 0 27 

SBS 1133+572 

40.76)(g 7 33 

41 17 +0 - 15 
^ a - a/ -0.15 

Mrk 1457 

41 74 + 0 - 96 
0.42 

4O.15*0 5 4 7 g 

2MASX J1 1570483+5249036 

40.67*g 

40.8 1*0 12 2 

2MASX J121 83945+4706275 

41 -68* q^2 

41 20 +0 - 33 
hi.zu_o3o 

2MASX J12384342+0927362 

41-69*0 42 

41 90 +015 
Hi.yu_o 12 

NGC 5695 

40.78*g 9 g 6 j 

9Q 7Q+0.42 


We plot L x , sf as a function of a x L[Nen] and L iA gn as a 
function of fi x L[Oiv] in Figure 3. We note that the abscissa 
errors are on the order of the symbol size, due to the 5(7 detection 
limit we imposed on IR detections (see LaMassa et al. 2012), and 
therefore are not plotted in Figure 3, though they are included in 
the errors for the Bayesian linear regression routine. Due to the 
large error on L v agn for NGC 7674, this source was dropped 
from the L^agn versus ft x L[o IV] analysis. 

In Figure 3, a dot-dashed line is plotted to indicate where 
the two quantities are equal while the gray shaded region 
illustrates the 3(7 confidence interval on the regression line. 
T x ,sf and L, A gn are approximately equivalent to a x L[Nen] 
and p x L|o iV |, respectively, albeit with wide scatter for a 
handful of individual sources. Deviations of the regression fit 
from equality at greater than the 3cr level are slight, and appear 
at higher luminosities for the star formation decomposition 
(>10 41 erg s 1 cm” 2 ) and moderate luminosities for the AGN 
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Table 9 

12 /xm Sample Star Formation and AGN Luminosities 


Galaxy 

log(G.SF) 

l°g(i X ,AGN) 

Comments 

NGC 0424 

40.99+° 0 ® 

40 72 +0 - 21 
hu./z_ 0 is 


NGC 1144 

40.78t°ri«i 

40.877° 0 2 2 \ 


NGC 1320 

40.33 + _‘ 0 ° 5 1 

40.427° 0 21 1 


NGC 1386 a 

72 Q 972+0.07 
■^•^—0.08 

39 S9+ 010 
j^.jz_oo7 


NGC 1667 

40.30+° 0 6 2 ° 1 

40.067°o 2 2 4 4 


F05 189-2524 

41 20 +0 - 78 
h-i.zu_ 030 

41 56 +015 

hi.^o_ 0 12 


NGC 3982 a 

7*0 49+0.10 

to -7 '7+0.04 
' —0.05 


NGC 4388 a 

40 52 +0 u 

< +u.jz_o i4 

39.877° 0 1 1 9 2 


NGC 450 l a 


39 40+0.09 
jv.h-u_o 09 


TOLOLO 1238-364 

40.63 + _° 0 7 5 2 7 

40 3 l +0 - 48 
hu.ji_ 045 


NGC 4968 

39 67 +0 ' 69 

jv.o/_0 33 

39.887° 0 1 1 2 2 


M-3-34-64 

41 47+0.45 
4 ' -0.21 

40 84 +0 - 60 

4+u.54f_ 060 


NGC 5 1 35 a 

40.867o.i9 

40.907° 0 ° 0 2 


NGC 5194 a 

72Q AQ+0.06 
JV.oy_Q 07 

39.087° 0 °o 4 4 


NGC 5347 

720 724+0.84 
— 0.48 

39.477° 0 1 1 5 2 


Mrk 463 a 

41 55 +0 - 41 
4fl.DJ_o 1 6 

41.697°o°o 7 9 


NGC 5506 

39.882°o 9 2 3 7 

41 -3870.09 

ObsIDs 0013140101, 0201830201, 0201830301, 0201830401 

NGC 5506 

0 

41 53+ 009 

ObsIDs 0013140201 and 0201830501 

NGC 5506 

40.067°o 4 2 5 1 

41.537o‘o6 

ObsIDs 0554170201 and 0554170101 

Arp 220 a 

40.512 0 o 1 o° 8 

39.587o°i4 


NGC 6890 

39.832° 0 f 9 

39.947°o 2 2 4 4 


IC 5063 a 

40 1 1 +009 
^ u - 11 -0.13 

40 69 +0 06 
HU.oy_0 06 


NGC 7130 a 

40.697o.o7 

40 94+ 0 - 05 
HU.yH-_ 004 


NGC 7172 

39.857o.36 

40.457° 0 °o 6 6 

ObsID 0414580101 

NGC 7172 

39.677°o 4 1 5 5 

40.187°o 1 1 2 2 

ObsIDs 0147920601 and 0202860101 

NGC 7582 a 

40 93+0-04 

4fU.ZJ_ 004 

39.727° 0 °o 3 4 


NGC 7674 

41 53+0.51 
■+1 -2 J_Q 21 

40.037 2 o 7 3° 9 



Note. a Sy2s observed with Chandra where the AGN was removed and extended emission was fitted. L x sv is derived from these flux values for these objects while 
L^agn is simply Lo.5-2keV _ £;t,SF- For the remaining sources, L x sf and L*,agn are derived as described in the text. 



Figure 2. Results of decomposing the soft X-ray flux (0.5-2 keV) into a star 
formation and AGN components, using the luminosity of the [Neil] line as 
a proxy for the former and the luminosity of the [O iv] line to parameterize 
the latter. The constants a and f5 were calculated using ordinary least-squares 
multiple linear regression, where we derive a = 1.07 and = 0.16. The 
overplotted line indicates where the two quantities are equal. Cyan triangles 
represent the 12 /xm Sy2s (with the X-ray variable sources connected by a 
vertical line), blue diamonds denote the subset of 12 /xm sources with high 
signal-to-noise Chandra imaging of the unresolved emission, and red squares 
mark the [O hi] selected sources. 

(A color version of this figure is available in the online journal.) 


parameterization (10 39 < L a A gn < 5 x 10 41 erg s -1 cm -2 ). As 
mentioned previously, the NLR for IC 5063 takes up a significant 
fraction of the X-ray extraction region, where contributions from 
resolved NLR emission could potentially bias the X-ray estimate 
of host galaxy star formation. However, as seen in Figure 3, 
where IC 5068 is the blue diamond at a x L[ NelI ] = 40.4 dex 
and L x sf = 40.1 dex, this source lies well within the shaded 
region. Hence NLR contamination appears to be negligible, 
possibly due to the removal of the AGN from the Chandra 
imaging analysis. 

We note that aperture bias can contribute to these slight 
disagreements: [Nen] ([Oiv]) is measured through the Short- 
High (Long-High) module on the Infrared Spectrograph on 
Spitzer, corresponding to sizes 4"7 x 11" 3 (ll'.'l x 22"3), 
whereas the X-ray apertures capture the full emission from 
the host galaxy (see Tables 1 and 2). Hence, biases can be 
introduced from comparing parameters that sample different 
physical scales of the host galaxy, introducing scatter into 
the relation. In the AGN parameterization, there are four 
obvious outliers from the relation: IC 0486 and Mrk 463 
(members of the [Om] sample with overpredicted L, A gn, 
red squares in Figure 3), NGC 4388 (from the 1 2 /xm sample 
with underpredicted L x ,agn: blue diamonds in Figure 3), and 
NGC 5506 (from the 12 /rm sample with overpredicted L x A qn, 
an X-ray variable source with a vertical line connecting the cyan 
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Figure 3. Left: L x sf vs. a x L[Nen]- Right: Z, Ai agn vs. /> x Zqoivp In both plots, the dash-dotted line indicates the line of equality, while the gray shaded regions 
illustrate the 3 cr confidence interval from Bayesian linear regression. Color coding is the same as in Figure 2. /. t sF and agn are derived from spectral fitting of the 
unresolved Chandra emission for the blue diamond data points while these parameters are estimated by scaling Lapec and L pow as noted in the text for the remaining 
data points. The soft X-ray and IR decomposition approximately agree at the 3 a level, with small deviations appearing at luminosities above 10 41 erg s -1 cm -2 for 
the star formation parameterization and between 10 39 < C.AGN £5 x 10 41 erg s _1 cm~ 2 for the AGN decomposition. 

(A color version of this figure is available in the online journal. 1 


triangles in Figure 3). Interestingly, there are no correlations 
between outliers in the L r agn parameterization and physical 
size of X-ray emitting region (where mismatch between aperture 
sizes between Chandra or XMM-Newton and Spitzer could 
introduce disagreements), quality of soft X-ray spectra, sample 
selection, source variability, or even method for estimating 
L t ,AGN- In general, however, deviations of the best-fit trends 
between L x versus a x i[Nen] and L.^agn versus /3 x Zqoiv] 
from equality at the 3er level are slight. Such an agreement 
demonstrates that our method of using low-resolution CCD 
spectroscopy to decompose the soft X-ray emission into a star- 
forming and an AGN component is moderately successful, with 
greater consistency in the star formation estimation. L r ,sF can 
then be used to estimate the host galaxy SFR in Sy2s using the 
Co. 5 — 2 kev calibrations from Ranalli et al. (2003) and Pereira- 
Santaella et al. (2011). 

5.1. Comparison with Far-infrared Derived SFRs 

For the 12 jim sample, we have tested the efficacy of our 
decomposition by deriving a soft X-ray SFR using the Pereira- 
Santaella et al. (2011) calibration and comparing this with the 
far-infrared (FIR) SFR using the Kennicutt (1998) calibration, 
using IRAS flux densities reported in Spinoglio & Malkan 
(1989). The [O m] sample, however, lacks a reliable independent 
host galaxy wide SFR. These sources were not detected by IRAS 
and the SDSS SFRs derived through the 3" fiber (see LaMassa 
et al. 2012 for a discussion) only covers the innermost regions of 
the galaxy, while the X-ray aperture captures the full emission. 
We therefore use only the 12 pm sample to test the X-ray- 
derived SFRs. 

To estimate SFR fir , we used the sum of 60 pm and 100 pm 
IRAS flux densities. We note that the Kennicutt (1998) calibra- 
tion is based on MIR through FIR luminosities for starburst 
galaxies, and in such systems, including the IRAS 12 pm and 
25 pm luminosities as part of the FIR luminosity is impor- 
tant. However, these MIR luminosities have significant contri- 
butions from dust heated by the AGN, so we only consider the 
60 pm and 100 pm luminosities when deriving SFRfir, which 
are less contaminated by the AGN. As shown in Figure 4, a good 
agreement exists between SFRo.5-2kev and SFR F ir, despite the 



log (SFR fir M 0 / yr) 


Figure 4. X-ray-derived SFRs (SFRo. 5 - 2 kev) using our estimate of L Xy sf and the 
calibration of Pereira-Santaella et al. (2011) as a function of FIR-derived SFRs 
(SFRfir), using IRAS 60 fim and 100 fim luminosities and Kennicutt’s (1998) 
calibration, for the 12 fim sample. The overplotted dashed line shows where 
the two SFRs are equal. A relatively good agreement exists, with the two most 
significant outliers with the largest SFRfir values (Arp 220 and F05 189— 2524) 
being ULIRGs. 

(A color version of this figure is available in the online journal.) 

latter SFR parameterization being most appropriate for starburst 
galaxies. This agreement suggests that any unresolved photoion- 
ized emission from the NLR does not introduce any systematic 
offsets into our results, but rather operates at a level less than the 
observed scatter. The two Sy2s that are the greatest outliers (i.e., 
largest SFR F ir values) are Arp 220 and F05 189— 2524, which 
are ultra luminous infrared galaxies (ULIRGs). 

6. CONCLUSION 

We have modeled the 0.5-10 keV spectra of two homoge- 
neous samples of Seyfert 2 galaxies to disentangle the starburst 
and AGN emission in soft X-rays (0.5-2 keV). Eleven of these 
Sy2s, observed with Chandra , had high signal-to-noise unre- 
solved emission after removing the AGN. We derive L a ,sf from 
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spectrally fitting this unresolved emission and assign L v A(;n to 
the difference between the total soft X-ray emission and L Xi sf. 

The remaining 24 sources were decomposed into AGN and 
star formation components by modeling the soft X-ray emission 
with a power law and thermal model. The luminosities of these 
sources were converted to L Xt sf and L a A gn values based on a 
scaling factor derived from XSINGS analysis on normal star- 
forming galaxies. To account for errors on the luminosities 
derived from spectral fits and the scaling factor, we executed 
Monte Carlo simulations, assuming a Gaussian distribution of 
random variables centered on the best- fit values of Lapec, T pow , 
and R, with the spread corresponding to the errors on these 
parameters. Using Equations (4) and (5), we calculate L v ,sf and 
L x ,agn from the 1000 simulated parameters and use the mode 
of this distribution as our estimate of the X-ray luminosities due 
to star formation and AGN activity. 

Our conclusions are summarized as follows. 

1. The soft X-ray spectra of 10 Sy2s were well fitted by two 
thermal model components, indicating the presence of a 
two-temperature gas. This result is similar to what has 
been observed in starburst galaxies (Dahlem et al. 1998; 
Strickland et al. 2004). 

2. For the subset of 1 1 Sy2s with high signal-to-noise Chanclrci 
imaging of unresolved host galaxy emission, estimates of 
the soft X-ray emission due to star formation and AGNs 
from both Chandra imaging analysis and from the Monte 
Carlo simulations agree at the 3er level. In the absence 
of resolved nuclear emission, scaling L pow and Tagn 
by the factors derived from XSINGS analysis is thus a 
reasonable method to estimate the AGN and star formation 
contributions to the soft X-ray emission. 

3. The independent decompositions of the soft X-ray lumi- 
nosity into star formation and AGN components using 
IR data as a proxy and scaling Tapec and L pow based 
on XSINGS galaxies largely agree within a 3a confi- 
dence interval. Deviations of the best-fit regression line 
from this agreement are slight and appear at higher lumi- 
nosities (>10 41 erg s 1 cm~ 2 ) in the star formation de- 
composition and at moderate luminosities ( 10 39 < L < 
5 x 10 41 erg s -1 cm -2 ) when describing AGN emission. 
Though more scatter in individual sources is evident in 
the AGN decomposition, the star formation relationship is 
more consistent among both methods. 

4. Comparison of our calculated X-ray SFR using the Pereira- 
Santaella et al. (201 1) calibration with an FIR-derived SFR 
from Kennicutt (1998) for the 12 /i m sample shows general 
agreement. 

We have demonstrated that analysis of low-resolution CCD 
X-ray spectra can effectively disentangle emission from AGN 
activity and star formation in 0.5-2 keV X-rays. Using the de- 
composition we have presented, L x ,sf can be used to cleanly es- 
timate the SFR in Sy2s using existing calibrations (e.g., Ranalli 
et al. 2003; Pereira-Santaella et al. 2011), complementing pre- 
vious studies in the optical and infrared and providing a more 
panchromatic view of the interplay between SMBH accretion 
and star formation. 
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